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Abstract 

Motivated by recent observations from Pamela, Fermi and H.E.S.S., we consider dark matter 
decays in the framework of supersymmetric SU(5) grand unification theories. An SU(5) singlet S 
is assumed to be the main component of dark matters, which decays into visible particles through 
dimension six operators suppressed by the grand unification scale. Under certain conditions, S 
decays dominantly into a pair of sleptons with universal coupling for all generations. Subsequently, 
electrons and positrons are produced from cascade decays of these sleptons. These cascade decay 
chains smooth the e"*" + e~ spectrum, which permit naturally a good fit to the Fermi LAT data. 
The observed positron fraction upturn by PAMELA can be reproduced simultaneously. We have 
also calculated diffuse gamma-ray spectra due to the excesses and compared them with the 
prehminary Fermi LAT data from 0.1 GeV to 10 GeV in the region 0° < / < 360°, 10° < |&| < 20°. 
The photon spectrum of energy above 100 GeV, mainly from final state radiations, may be checked 
in the near future. 
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I. INTRODUCTION 



Electron, proton, photon, neutrino and their antiparticles are stable, at least on the 
cosmological time scale. Detection of these particles from cosmic rays provides an interesting 
window to look into the deep universe. Recently, the PAMELA experiment reported a 
significant excess in the positron fraction e'^ /{e^ + e~) between 10 GeV and 100 GeV [1]. 
On the other hand, the measured antiproton to proton flux ratio appears to be consistent 
with predictions [2]. More recently, the Fermi LAT collaboration observed a smooth e"*" + e~ 
spectrum with high accuracy. It is found to be falling as E~^'^ from 20 GeV to 1 TeV 
[3], much harder than the predictions of conventional models. The H.E.S.S. collaboration 
measured the e"*" + e~ spectrum from 600 GeV up to several TeV [4] , which is consistent with 
the Fermi data in overlapping regions and steepens at about 1 TeV towards higher energy. 

These excesses of electrons and positrons could be due to unidentified astrophysical 
sources, e.g., nearby pulsars or supernova remnants [5-7]. However, an explanation via 
dark matter (DM) annihilation or decay is, arguably, a much more interesting possibility, at 
least from the perspective of particle physics. The electron and positron spectra alone, even 
with higher precision and broader energy range, cannot decisively decide which explanation 
is more plausible [8]. Hopefully, the energy spectrum and the angular dependence of cosmic 
gamma rays [9-13], to be measured by the Fermi LAT in the near future, may provide a 
more definite answer. For the DM interpretation, the mass of the DM should be around 
several TeV, to provide the excesses from 20 GeV to 1 TeV and steepen sharply above 
1 TeV. Furthermore, traditional WIMP DM candidates usually produce extra antiprotons. 
As Pamela does not observe any deviation on antiproton spectrum from the anticipation, 
WIMP DMs are now disfavored as potential sources of the observed cosmic-ray excesses. 
Still, there are plenty of freedoms for both DM annihilation and decay to reproduce the 
experimental spectra reasonably [14, 15]. For DM annihilation, a large boost factor in 
the order of 10^ to 10'^ is needed for the theory to be consistent with the relic abundance 
measured by the WMAP [16]. As the dumpiness property of the DM distribution falls far 
short of such a large factor, one usually resorts to nonperturbative Sommerfeld [17-20] or 
Breit-Wigner [21-23] enhancement in model buildings. For DM decays, the lifetime should 
typically be around the order of lO^^s to fit the data[14, 24, 25], which is much longer 
than the lifetime of the universe. Therefore the DM decay rates will not affect the relic 
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abundance appreciably. 

The energetic flux produced from DM annihilations/decays would inevitably emit 
gamma rays. These gamma rays depend on the DM density as for annihilations and p for 
decays. This will lead to different angular dependence of the gamma ray spectrum, which 
may be measurable in the near future to differentiate these two scenarios. The gamma ray 
spectrum can also be used to differentiate DM explanations from astrophysical ones. 

In this paper, we will focused on DM decays. Notice that a lot of suppression will be 
needed for a TeV scale particle to have a lifetime ~ lO^^s. If it decays via dimension 
four operators, tremendous fine tunings will be needed. If it decays via dimension six 
operators, it still needs to be suppressed by a scale ~ 10^^ GeV, which turns out to coincide 
with the grand unification theory (GUT) scale [26-28]. In the same spirit of Refs. [29, 
30], we will take a singlet as the dark matter candidate and provide a detailed analysis 
in the frame of supersymmetric SU(5) GUT. To be consistent with the Pamela antiproton 
measurement, squark masses are assumed to be heavier than that of the SU(5) singlet S, 
so the S decay would be quark phobic. The S then decays dominantly into slepton pairs 
with a universal coupling for all generations. These sleptons decay quickly into leptons and 
lightest supersymmetric particles (LSPs), if R-parity is conserved. In this framework, we 
have obtained a reasonable fit to all data from Pamela, Fermi and H.E.S.S.. 

The fluxes from S decays are inevitably accompanied by hard photons: coming from 
final state radiations (FSR) of cascade decays, including S*— t-t— t-t— t-tt^— 7-27 and the 
inverse Compton scattering (ICS) on the interstellar radiation field (ISRF). The gamma ray 
fluxes could have Galactic and extragalactic origins. We have calculated all these gamma 
ray spectra and compared them with the recent Fermi LAT measurement in the region 
0° < / < 360°, 10° < |6| < 20° [31]. 

This paper is organized as follows. The supersymmetric SU(5) model plus a singlet S 
is presented in Section II, where we have also discussed the possible decay channels of S 
in some detail. In Section III, a reasonable fit is obtained to reproduce the observed 
fluxes, by tuning relevant parameters in the model. Section IV is devoted to the study of 
gamma-ray spectra from excesses. Finally we conclude with a summary in section V. The 
component field structure of the dimension six effective operators will be presented in the 
Appendix. In this paper, we have used the Navarro-Frenk- White (NFW) halo model [32] for 
DM distribution and the MED propagation model [33, 34]. For other halo and propagation 
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models, the conclusions are similar. In addition, all computations on astrophysical effects 
are performed semi-analytically instead of using the GALPROP program^. 

II. A SUPERSYMMETRIC SU(5) MODEL 

If the observed excesses come from DM decays, the lifetime of a TeV scale DM should be 
~ lO^^s. Such a long lifetime can be naturally realized through decays via GUT suppressed 
dimension six effective operators, similar to proton decays. This provides a strong motivation 
to study DM decays in the framework of grand unification theory [29, 30, 35-38]. 

In the minimal supersymmetric SU(5) model, the dark matter candidate would be the 
LSP, which is absolutely stable if R-parity is conserved. In addition, the mass of LSP is nor- 
mally around several hundred GeV, which is too small to account for the Fermi and H.E.S.S. 
data even if it decays. To make a minimal extension, one can introduce an SU(5) singlet 
S as the dark matter candidate^ [29]. As S is neutral in the standard model (SM) gauge 
group, it does not disturb the gauge coupling unification. To eliminate lower dimensional 
operators which may lead S to decay too fast, we impose a Z2 symmetry on the theory, 
under which S is odd while all other particles are even. Then 5* can decay into the MSSM 
particles only through dimension six operators, suppressed by MQ^rp. Assuming R-parity 
conservation and the Z2 symmetry, all possible dimension six operators are [29]: 

^GUT ^GUT ^GUT ^GUT 

Some or all of these operators may appear at the TeV scale when one integrates out heavy 
particles of the GUT scale. Here Wa are the supersymmetric field strengths of SM gauge 
groups, Hu and Hd are the chiral superfields for Higgs, 5 and 10 are anti- fundamental and 
antisymmetric tensor representations of SU(5), respectively 

5^ = (rf^d^c^^e,-^.)i (2) 

^ Web page: http://galprop.stanford.edu/web_galprop/galprop_home.html 

^ If R parity is conserved, the neutralino LSP would also be part of the DM. But for simplicity, we assume 
here that S is the dominant component of DM and the LSP contributes just a small portion to the relic 
density. We will show that such a scenario is feasible in the next section. 
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Operators S'^ SH^f^^~^Hu{d) and S'^ SWaW^ in Eq.(l) may lead to final states containing 
significant number of quarks or mono-energetic gamma ray lines. To be consistent with 
experimental observations, these operators should be further suppressed, most likely due to 
unknown physics at the GUT scale. They will simply be neglected from now on. Operators 
S+STl and 5+5Tr(10+10) in Eq.(l) may be rewritten in the form 



/l/f2 



(4) 



Here the summation is over all lepton and quark chiral superfields. Assuming the singlet 
scalar develops a vacuum expectation value (VEV) < s' >, the Z2 symmetry is spontaneously 
broken and both components s) in S will decay. 

Expanding Eq. (4) in terms of component fields, one has 

— < s > T{d^ipa^ip) + i <s> i)*{d^ipa^s)+ <s> s*^*D^^ + h.c. + ... . (5) 

Here we have dropped total divergence terms. Operators from F-terms have also been 
neglected as they are suppressed by the leptonic Yukawa coupling constant. In addition, 
these operators will lead to many body decays which are further suppressed by phase spaces. 
Details of the expansion will be provided in Appendix A. 

To fit the e+ fluxes data which steepens sharply above TeV, the mass of s and 's will be 
assumed to be around several TeV. In addition, the squark masses are assumed to be heavier 
than the DM mass while the slepton masses to be about several hundred GeV. So, s and 
's can only decay into leptons, quarks and sleptons, and have no squarks in the final state. 
The assumptions on squark and slepton masses seems to be plausible, because squarks are 
much heavier than sleptons in general. This is due to the fact that squarks are color charged, 
which may affect drastically the renormalization group equations for the squark masses. 

The decay width of s (s) due to the first two operators in Eq.(5) is proportional to the final 
state quark or lepton mass square. For the first operator 's*{d^ipa'^ilj), the dominant decay 



5 



channel is J—)- ti, which is suppressed by M^/M^ compared to the third operator 's*ip*\I\ip 
in Eq.(5). The main decay channel of s through the second operator ■ip*{dfj_ipa'^s) is s — rr, 
which is again suppressed by M^/M^ compared to the third operator in Eq.(5). Thus the 
decays of s will not be considered. The DM 's decays dominantly into a pair of sleptons, 
with universal coupling for all generations. For simplification, we will simply neglect the 
first two operators in Eq.(5) and only consider the operator s'*?/'*n?/' in the following. The 
remaining operator can be further rewritten as 



J2t7^ <s>r{lLDlL + iRDln) (6) 
with I = e, Jl and r. The corresponding decay width reads 



T~ = -V I - (7) 

Taking Mqut = 10^^ GeV, M, ~< ? >~ a few TeV and Mj ~ several hundred GeV, the 
lifetime of 's would be around lO^^s, as one has hoped. Notice also that the decay width is 
proportional to M~, so that even slightly different masses between e, Jl and r may lead to 
very different branching ratios. 

With R-parity conservation, the slepton would decay to the LSP and lepton quickly.^ 
can be produced through the following cascade decay chains: 

• selectron chain: S' — > e' — )■ e 

• smuon chain: 's ^ Jl ^ fi ^ e 

• stau chain: s — r — r — ;> e 

In total, the fluxes due to DM decays at the source are 

Q?"ir,E) = ^ ' ^J > J^' (8) 

I 

Here the summation is over all three cascade decay chains. T-^^ is the decay width of the 
/ cascade decay chain and M^*^ is the DM mass. Since the lifetimes of sleptons, muon 



^ Here the LSP is assumed to be the neutrahno. If the LSP is the gravitino and s is heavier than s, s would 
decay into s predominately, instead of the sleptons. 
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and tau are extremely short compared with the DM decay, we can take the approximation 
r~"^^ = Tj. And dN~^^ /dE is the spectrum of electron or positron per DM decay via 
a particular / chain. For the stau chain, the e"^ spectra are obtained by using PYTHIA 
package [39]. p^^\r) is the DM mass density which is model-dependent. As an illustration 
we adopt the NFW halo model [32] 

r V 1 + "^rs J 

with solar system position Tq = 8.5 kpc, the DM density at earth = 0.3 GeV/cm^ and 
Ts = 20 kpc. 



III. ELECTRON AND POSITRON EXCESSES FROM DARK MATTER DECAY 
A. Positron and Electron Propagation 

Shown in Eq.(8) are the fluxes due to DM decays at the source. However, only the 
e"^ fluxes at the Earth are observable. It is thus necessary to consider the propagation of 
electrons and positrons in the Galaxy. The flux per unit energy at an arbitrary space-time 
point is given by 

$r(t,f,E) = ^/f^^(t,f,E). (10) 

For energetic e^'s that we consider, their velocity is approximately equal to the light 
speed c. The function f^^'\t,r,E) satisfies the diffusion-loss equation 

^ = K{E) . V'lr + ^ (5(i?)/e^^) + Qr • (11) 

Here the convection and advection terms have been neg lected. Qf^^ is due to the DM decays 
as given in Eq.(8). K{E) stands for the diffusion coefficient which is related to the rigidity 
of the particle. For e^, it can be parameterized as 

K(E) = KoiE/GeVf . (12) 

B{E) = E"^/ (GeV-TE) is the effective energy loss coefficient with te = lO^^s, which describes 
the energy loss of due to ICS on the ISRF and synchrotron radiation. Eq. (11) can be 
solved in a solid fiat cylinder parameterized by {r,z,6), with z G [—L,L] in z direction 
and r G [0, 20 kpc] in radius. The solar system corresponds to the position {tq, Zq,9q) = 
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0.806 


0.891 


0.721 
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0.071 



Table I: Parameters in Eq.(15) in accord with the NFW halo model and MED propagation model. 

(8.5 kpc, 0, 0). /(P^(t, X, E) is assumed to vanish on the surface of the fiat cyhnder, which 
serves as the boundary condition for this equation. For the MED propagation model, L, S 
and Kq are chosen to be 4 kpc, 0.70 and 0.0112 kpc^/Myr, respectively [40]. 

If f^'^'\t,x,E) does not change with time, a steady state solution Eq. (11) can be 
obtained semi-analytically [34, 40-42]. The fluxes at the Earth are 
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Here X£){E,E') describes the diffusion length from energy E' to E, which can be parame- 
terized as 



^D = 4Aori5 I 



The function /(A/j) is given by: 

H^d) = tto + ai tanh 
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(15) 



with / = log]^Q(Az)/kpc). I{Xd) contains the whole information of the NFW halo model and 
MED propagation model. Parameters in Eq.(15) have been estimated numerically in Ref. 
[40] and listed in Table I. 

B. Positron and Electron Backgrounds 



For interstellar background fluxes of e^, we use the "model 0" proposed by the Fermi 
LAT collaboration[7], which can be parameterized as [43]. 

82.0e-o-28 



K-'{E) 



1 + 0.224e2-93 ' 
38.4e-^-^8 

1 + 0.0002e5-63 



24.0e 



-3.41 



(16) 
(17) 



in units of GeV m ^s ^sr ^, where e = E/lGeV. 
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For electron/positron fluxes at the top of the Earth's atmosphere $®t, solar modulation 
effects should be considered. Adopting the force fleld approximation, one has [44] 

^^E^) = ^mEis) , (18) 
^is 

where stand for interstellar fluxes and = Ejs + \Ze\(j)F, with (pp = 0.55GV as a 
typical value. It is clear that, at energies larger than lOGeV, solar modulation effects could 
be neglected as E^ ^ Ejs- 

Finally, the fluxes and positron fraction at the top of the Earth's atmosphere could 
be expressed as 

El X <|.e(i?e) =Elx + $,^(^,5) + N ■ ^'Jl'iEjs) + ^f-'\Ejs)] , (19) 

Eis 

^e^iE^) K''iEis) + ^f.''{Eis) ^20) 



$e+(i5e) + *e-(i?e) ^l^'iEjs) + ^f^'iEjs) + N ■ ^I'JiEjs) + ^f^''{Ejs) ' 
Here is a normalization factor standing for the uncertainty of the electron flux. In this 
paper = 0.8 is chosen to flt the experimental data. 



C. A Fit of PAMELA and Fermi LAT Data 

For illustration, we choose DM mass Mj = 6.5 TeV, Mqut = 10^^ GeV, < s >= 20 
TeV, Me = 380 GeV, = 370 GeV, = 330 GeV and Mlsp = 300 GeV.^ With 
this parameter set, the decaying DM produces extra fluxes from 10 GeV to 1 TeV, as 
shown in Fig. la. The main contribution comes from the selectron chain. The cascade decay 
? — )■ e"^ + e~ —^e^ + e~ smoothes the + e" spectrum and naturally allows for a good flt to 
the Fermi LAT measurement. The fluxes steepen above ITeV sharply, which is consistent 
with the H.E.S.S. observation. The positron fraction are shown in Fig.lb, compared to the 
data of PAMELA, AMS-01, CAPRICE and HEAT. 

^ We have checked expUcitly that, for this set of sleptoii masses, the neutrohiio LSP could be only a minor 
part of the DM. For example, by using DarkSUSY package[45], we have obtained ^l^sph^ = 0.009, with 
the neutralino mass spectrum Mlsp = M^i = 300 GeV, M^^ = 315 GeV, M^^ = 630 GeV, M^^ = 690 
GeV and the gaugino fraction to be 0.03. 
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Figure 1: The SU(5) model gives a reasonable fit to the Fermi, H.E.S.S. and Pamela data with the 
example set of parameters given in the text. Left: Decaying DM produces extra e~ + e"*" fluxes 
above background via three different cascade decay chains. The green dot line, black dash dot 
line and pink dash dot dot line represent the selectron chain, the smuon chain and the stau chain, 
respectively. The red solid line includes all contributions from our fit. The black dash line shows 
the background as discussed in the text. Right: Including the e"*" background flux, decaying DM 
predicts a positron fraction which fits the experimental data. The red solid line shows the result of 
our fit, and the black dash line indicates the background. 

IV. DIFFUSE GAMMA-RAYS FROM THE EXCESSES 

The excesses are inevitably accompanied by photons coming from the FSR, ICS and 
synchrotron radiation stemming from them. 

(1) FSR: The bremsstrahlung of fluxes leads to the emission of energetic photon flux 
^FSR- Moreover, in our model the stau chain contains r lepton which emits hard photons 
via the process r — t- tt*^ — 7 + 7. This mechanism is significant, especially at the high 
energy end of the spectrum. The largest energy of FSR photons could be around 

which could be probed by the H.E.S.S. collaboration. Notice also that the spectrum of FSR 
is quite model-dependent. In addition, the FSR could come from within or without our 
Galaxy (<f^^y$|f«). 

(2) ICS: The ICS radiation is produced when the excesses scatter on the ISRF. In 
our Galaxy the ISRF includes the cosmic microwave background (CMB), star light and 
the infrared light, while outside of our Galaxy the CMB component is dominant. The 
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corresponding ICS fluxes ^fcs ^^"^ ^fcs^ which should be observable by Fermi LAT, are 
closely related to the excesses. Therefore the ICS spectrum is, to some extent, model 
independent as long as the DM model can reproduce the Fermi and H.E.S.S. spectrum 
with reasonable accuracies. 

(3) Synchrotron radiation: During propagation, the fluxes radiate photons in the 
Galactic magnetic fields. These photons should be very soft (around 10~^eV) and outside 
the energy range explored by Fermi LAT and H.E.S.S. experiments. They will not be 
considered in the following. 

Notice that the extragalactic gamma rays are roughly of the same order as the Galactic 
ones. But the extragalactic component is isotropic while the Galactic one has angular 
dependence. The total gamma ray flux is obtained by summing all these contributions: 

% = + + $f + ^fSs (21) 

Specifically, we consider only photons in the region 0° < / < 360°, 10° < \h\ < 20° in the 
following, as Fermi LAT has released the data in this region recently [31]. 



A. Galactic Gamma Rays from FSR 

As photons propagate almost freely in the Galaxy, the differential flux of photons received 
at the Earth in a given solid angle dQ is given by [40] 
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P0 jsr^Y (22) 



dE^dn AttM^^ ^ ' dE^ 



Here the factor of 2 takes into account the fact that both leptons and anti-leptons contribute 
equally to the FSR flux of gamma rays. dN~^'^ /dE^ is the photon spectrum per DM decay 
via a specific slepton chain. PYTHIA package [39] has been used here to obtain these 
spectra. J encodes all the astrophysical information which is defined as 

jAn= [dn r^P^, (23) 

J Jo rQ Pq 

where the parameter s is integrated along a line of sight. The parameter s can be linked to 
parameters r, I and b by 

r (s, /, &) = Y^s^ + r| — 2srQ cos / cos b (24) 
11 
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Figure 2: Galactic gamma ray spectra from FSR (left) and ICS (right) via three different decay 
chains, plotted in the region 0° < / < 360°, 10° < |fe| < 20°. 

From Eq.(23), one obtains 7 = 2.4 in the region 0° < / < 360°, 10° < \b\ < 20°. The 
photon spectra from Galactic FSR are plotted in Fig. 2a, which peak around several hundred 
GeV. Notice that the stau chain gives a large contribution to the photon spectrum due to 
r — ;> 7r° — 7 + 7. 



B. Galactic Gamma Rays from ICS 

A pedagogical review about ICS was provided in [46]. We will calculate the ICS gamma 
rays semi-analytically, following Refs. [12, 43, 47, 48]. 

The differential flux of ICS photons received at the Earth in a given solid angle dQ with 
energy between and + dE-y can be expressed as: 

^(hGAL 2r r r°° r°° /•A/^'^Va dn^^^iE A 



(26) 

here /e(r, E^,) denotes initial electron number density and fisRpi'f^i the ISRF photon num- 
ber density. The factor of 2 reflects that both electrons and positrons contribute to the ICS 
gamma rays equally. The Compton cross section is given by the Klein-Nishina formula 

' 2glng + l + g-2g +— — — — (1 -g) (26) 



dE^ A-ile ^ ^ ^ ■ ■ ^ " ' 2(1 + gF) 

where 

E^ 47ee _ Ee 



r = 7e = — . (27) 



T[Ee — E^) ' me ' m, 
12 



Here (Tt = 0.67 barn is the Compton scattering cross section in the Thomson hmit and rrie 
is the electron mass. Kinematics requires that e < E'^ < (l/E'e + 1/47^6)"^. 

The initial electron or positron number density /e(r, E'g) can be obtained by solving 
Eq.(ll) at each position. Notice that Eq.(ll) is dominated by the energy loss term at high 
energy. That is to say, can not propagate far from the production position before losing 
most of their energy. Therefore, Eq.(ll) may be solved point by point approximately 

i''-(''^') = ^)W^T.^nE.) (28) 

with 



To approximate further, we will assume that the ISRF photons have the same energy 
spectra at any point in the region 0° < / < 360°, 10° < \b\ < 20°. That is to say, the 
number density of ISRF fisRpi'f^, = fiSRpi^), which can be described by three blackbody- 
like spectra roughly [47]: 

fiSRFie) = Y: ^^^27^1 (30) 

j=l,2,3 

with Ti = 2.753 K, iVi = 1 for the CMB, T2 = 3.5 x lO'^ eV, N2 = 1.3 x 10~^ for the 
infrared light and Tg = 0.3 eV, N3 = 8.9 x lO^^^ for the star light. 

In order to separate astrophysics and particle physics information, Eq.(25) can be rewrit- 
ten as 



dEydil 



^'J^J^L^l ^^^^ ^^^^E^^fryjiEe)fisnAe). (31) 



Shown in Fig.2b is the ICS photon spectra in our model. One sees that the gamma ray 
fluxes come mostly from the selectron chain and steepen sharply above 1 TeV. 



C. Extragalactic Gamma Rays from FSR 



To study gamma rays from the outside of our Galaxy, the effects due to the expansion 
of the Universe should be considered. By turning the line-of-sight integral into a redshift 
integral, the differential flux of isotropic photons of the extragalactic origin is given by 
[12, 49] 



dE^dVl 



An Mdm ^ 
I 



Err 



Hiz) 



(32) 
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Figure 3: Extragalactic gamma ray spectra from FSR (left) and ICS (right) via three different 
decay chains, plotted in the region 0° < / < 360°, 10° < |6| < 20°. 

Here z is the redshift, li{z) = //qa/^a + ^m{z + 1)^ is the Hubble expansion rate, with 
Hq = lOOh kms~^Mpc~^ and the present day normahzed Hubble expansion rate h = 0.72[16]. 
Pc = 5.5 X 10~^GeV/cm'^ is the critical density of the Universe. We also take the dark matter 
density Qdai = 0.21, the dark energy density = 0.74 and the matter density Qm = 0.26 
[16]. The spectrum dN~^^ /dE'^ is the same as that in Eq (22), except that the redshift 
effect has been included. The parametric form for the optical depth t{E^,z) of the "fast 
evolution" model could be found in Refs. [50, 51]. Fig. 3a shows those contributions of our 
model. Again the stau chain is important here because of the 7r° channel. 



D. Extragalactic Gamma Rays from ICS 

We adopt a semi-analytical calculation following Refs. [48, 49]. Concerning the dilution 
effect due to the expansion of the Universe, the diffusion-loss equation of electrons and 
positrons becomes 

a fDM a fDM a 

\- = H{z)Ef-^ + ^ [B^^iz^ E^)r^^] + Qf ^ . (33) 

Here the extragalactic energy loss rate B^^{z^ E) is given as 

5^^(^, Ee) = ^aT-ilpcMB{l + zf (34) 

with pcMB = ^■yPc = 0.26 X 10~^GeV/cm^ the present-day CMB energy density. For 
energy around several hundred GeV, the timescale of energy-loss is E / B^^[z, E) ~ lO^^s, 
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which is much less than the Hubble time. That is to say, basically do not feel the 
redshift effect before losing most of their energy. So the Hubble term in Eq(33) can be 
safely neglected. The e"^ spectrum from DM decay can then be solved as 

^°"(^- ^''> = J^S^^Y.^r'yiiE.) ■ (35) 

Finally, the differential flux of extragalactic ICS photons received at the Earth in an arbitrary 
solid angle dQ with energy between and + dE^ can be expressed as: 



d^f^s _2c dzdedEe j^^r . da^^^{z, E,,e) 

Jo JlVzfWr^'''^^' ^"'^ 



dE^dn in Jo {1 + zfHiz)' " ' ' dE'^ 



(36) 



The spectrum f-y{z, e) of the background CMB radiation at redshift z is given as 

1 



with T = 2.753 K. The photon spectra from extragalactic ICS are plotted in Fig. 3b, which 
are dominated by the selectron chain contribution. Because of the redshift, the spectrum 
drops rapidly at high energy. 

Finally, Fig. 4 shows the total gamma ray spectra including all contributions. One can 
see that the FSR 7-rays dominate at higher energies while the ICS ones dominate at lower 
energies. Extragalactic gamma rays are not significant at high energy due to the redshift 
effect. The total gamma ray spectrum from excesses are consistent with the preliminary 
Fermi LAT data [31] from 0.1 GeV to 10 GeV, as shown in Fig. 4. The predicted gamma 
ray flux around several hundred GeV may be tested by the Fermi satellite in the near future. 



V. SUMMARY 

In this paper we have studied the DM decay in supersymmetric SU(5) models. An SU(5) 
singlet S, instead of LSP, is assumed to be the dominant component of DM. With R-parity 
conservation and a spontaneously broken Z2 symmetry, the singlet S can decay into visible 
particles through dimension six effective operators suppressed by the GUT scale. Assuming 
the squarks to be heavier than S, S decays dominantly into a pair of sleptons through 
the effective operator 's*{Il DIl + ^rOIr ). Typically, the lifetime of S is around 10^^ s, 
much longer than the age of the Universe. Since the decay products of S do not contain any 
quarks, our model is consistent with the Pamela antiproton measurement automatically. For 
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Figure 4: All gamma ray spectra in the region 0° < / < 360°, 10° < |6| < 20°. Dots stand for the 
preliminary Fermi LAT data [31]. 

illustration, we have chosen M^^ = 6.5 TeV, Mqut = 10^^ GeV, = 380 GeV, = 370 
GeV, = 330 GeV and M^sp = 300 GeV. With this parameter set, we have shown that the 
S decays can account for the PAMELA, H.E.S.S. and Fermi LAT excesses. Numerically, 
we have adopted the NFW profile for the dark matter distribution and the MED propagation 
model for the cosmic ray propagation. To interpret these results, one should keep in mind 
that there exist substantial astrophysical uncertainties about background, propagation 
and the DM distribution. 

When propagate from the decay position to the Earth, hard photons are emitted 
inevitable due to inverse Compton scattering and final state radiation. Future measurements 
of the diffuse gamma ray may distinguish DM explanations from astrophysical explanations 
by looking into the energy and angular distributions. We have calculated the gamma ray 
spectra in our model. The predicted photon spectrum are compared to the preliminary Fermi 
LAT measurement from 0.1 GeV to 10 GeV in the region 0° < / < 360°, 10° < \b\ < 20°, 
which seems to be consistent with each other. The total gamma ray spectrum are dominated 
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by photons from Galactic final state radiation for the photon energy above 100 GeV, which 
may be tested by Fermi LAT in the near future. 
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Appendix A: the component field structure of S'^S'<I>^<I> 

We now provide the component field structure of the dimension six operator S'"*'5'$^<1>. 
Define 



Hy) = Hy) + V2ei^iy) + o'fm, $+(2/+) = r{y+) + V2dtp{y+) + e^F;{y+) (ai) 

Siy) = ~siy) + V2esiy) + e'F.iy), = + V2^s(y+) + e'F:iy+) (A2) 

with = + zOa^'e and ?/+™ = x"' - zOa'^e. 



Products of chiral superfields are again chiral superfields, and likewise for their conjugates. 



Define again 



Aiy) = ^y)Siy) = ~a{y) + V29a{y) + O^FM 



(A3) 



(A4) 



The component fields of these composite superfields are 



a = tjjs, a = tjjs + s4>, Fa = 4>Fs + sF^ — ijjs 



(A5) 



(A6) 



The 6669 term of the dimension six operator 



($5)+ ($5) 




(A7) 
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is given as [52]: 

F*Fa + ^rna + hla*^ - ]^d„Xd'^a + ^-draaa'^a - ^aa^'d^a (A8) 

Assuming the singlet scalar develops a vacuum expectation value (VEV) < s' >, the Z2 
symmetry is spontaneously broken and both components (s, s) in S will decay. Eq. (AS) 
can be reexpressed via Eq.(A5,A6). 

^ (i<I> s*{d^i)a^^) + i<s> ^*{d^i)(jn)+ <s> T^*U^) + h.c. + ... (A9) 

Here we have dropped total divergence terms. Operators from F-terms have also been 
neglected as they are suppressed by the leptonic Yukawa coupling constant. In addition, 
these operators will lead to many body decays which are further suppressed by phase spaces. 
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